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Topological excitations are found throughout nature, in proteins and DNA, as dislocations in 
crystals, as vortices and solitons in superfluids and superconductors, and generally in the wake of 
symmetry-breaking phase transitions 1 . In fermionic systems, topological defects may provide 
bound states for fermions that often play a crucial role for the system's transport properties. Fa- 
mous examples are Andreev bound states inside vortex cores [2], fractionally charged solitons in 
relativistic quantum field theory [3| H] , and the spinless charged solitons responsible for the high 
conductivity of polymers 5 . However, the free motion of topological defects in electronic systems 
is hindered by pinning at impurities O |6]. Here we create long-lived solitons in a strongly in- 
teracting fermionic superfluid by imprinting a phase step [THS] into the superfluid wavefunction, 
and directly observe their oscillatory motion in the trapped superfluid. As the interactions are 
tuned from the regime of Bose-Einstein condensation (EEC) of tightly bound molecules towards 
the Bardeen-Cooper-Schrieffer (BCS) limit of long-range Cooper pairs 10 , the effective mass of the 
solitons increases dramatically to more than 200 times their bare mass. This signals their filling with 
Andreev states pTHTSj and strong quantum fluctuations 14-19 . For the unitary Fermi gas, the 
mass enhancement is more than fifty times larger than expectations from mean-field Bogoliubov-de 
Gennes theory [TJl [20] . Our work paves the way towards the experimental study and control of 
Andreev bound states in ultracold atomic gases [111 121] . In the presence of spin imbalance, the 
solitons created here represent one limit of the long sought-after Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state of mobile Cooper pairs [2X1122] . 



The emergence of effective particles from an existing 
background of matter waves is one of the key concepts 
in physics, from Landau's quasi-particles to the fraction- 
alised charges in quantum Hall states, the Higgs mode 
in elementary particle physics and cosmic strings in the 
early universe [P. Superfluids are a paradigmatic form 
of such quantum matter. They are described by a com- 
plex macroscopic wavefunction that is rigid against twists 
of its phase. The ground state of the superfluid thus 
has uniform phase, and small perturbations propagate 
as sound waves. A slower, non-linear excitation - the 
dark soliton - occurs when the phase is twisted substan- 
tially over a short range. In the extreme case of a phase 
jump by 180 degrees, the wavefunction changes sign and 
crosses zero at the location of the jump, creating a sta- 
tionary black soliton. In weakly interacting BECs all 
bosons reside in the condensate, so the particle density 
vanishes at a black soliton, and is reduced for a moving 
dark soliton. Dark solitons in BECs have been created \ Ty 
[9l [23H26] and studied extensively theoretically [18 . They 
are well described by solutions to the non-linear Gross- 
Pitaevskii equation. In fermionic superfluids, solitons are 
phase twists in the wavefunction of fermion pairs [TTl [27] . 
In the limit of tight molecular pairing, the molecular con- 
densate is still described by the Gross-Pitaevskii equa- 
tion and stationary solitons are again devoid of particles. 
However, in the BCS regime of long-range overlapping 
Cooper pairs, only a minute fraction of particles near the 
Fermi surface takes part in pairing, and the reduction of 
the pair wavefunction at the soliton affects the density 
only very weakly. Solitons in BCS superfluids should 



thus be filled with normal fluid. Indeed, the depletion of 
the pair wavefunction gives rise to Andreev bound states 
localised at the soliton (TTJ [22] . In three dimensions the 
soliton is a planar defect, one therefore expects a two- 
dimensional Fermi gas of Andreev states to be confined 
to the soliton [TTl [22] . 

In the limit of weak attractive interactions, there exists 
a direct connection between solitons in a BCS superfluid 
and those studied in relativistic quantum field theory [3] 
and conducting polymers [5 . Here the size of Andreev 
bound states becomes much larger than the interparticle 
spacing, and their energy much smaller than the pairing 
gap. In this limit, the Bogoliubov-de Gennes equation for 
the bound state simplifies to the Andreev equation [22] , 
a Dirac equation where the pair wavefunction A(z) plays 
the role of a spatially varying mass coupling particles and 
holes. For vanishing transverse momentum it reads: 

Here, cFx.z are Pauli matrices, vp = is the Fermi 
velocity, kp = (37r^n)^/^ is the Fermi wavevector re- 
lated to the total density n, UnC^^^^ and v^e'^^^^ are 
the quasi-particle amplitudes for particles and holes and 
En is the quasi-particle energy. When A (2:) changes 
sign, this equation allows for one normalisable zero- 
energy mode per spin state, localised at the zero crossing 
of the gap, which famously carries half-integer particle 
number [3 and gives rise to the conductivity of poly- 
mers [5]. A self-consistent solution for the pairing field 
is [5] A(z) = Aq tanh(2:/^), where Aq is the bulk pairing 



gap and ^ = is the BCS coherence length. This soh- 
ton solution is shown in Fig. la), along with the density 
of atoms in the bound state, demonstrating its localised 
nature. In a BCS superfluid, the Andreev bound state en- 
ergy is never strictly zero owing to the non- vanishing ki- 
netic energy cost /m(^^ ~ A^/Ep of forming the bound 

state [11122] 5 where Ef = ^^^^ is the Fermi energy. 

If A{z) was restricted to be real, as in the theory of 
conducting polymers [5 , solitons would be topologically 
protected. However, at finite temperature solitons can 
accelerate due to dissipation, and the pairing field ac- 
quires an imaginary part. When the soliton reaches a 
critical velocity, it is expected to decay into sound waves 
or pair excitations [28] . Apart from this thermody- 
namic instability, in three dimensional systems the soli- 
ton can decay into vortices via the so-called snake in- 
stability |91i23J. Finally, for strong interactions between 
fermions, it is a priori not obvious that solitons are stable 
against quantum fiuctuations p^HTO] . 

Here we create and observe long-lived solitons in a 
strongly interacting fermionic superfiuid of ^Li atoms 
near a Feshbach resonance. Solitons are created via phase 
imprinting (see Fig. lb), a technique successfully em- 
ployed for weakly interacting Bose condensates [71 [3 [25] . 
The superfiuid containing typically ~ 2 x 10^ atom pairs 
is prepared in an elongated trap with cylindrical symme- 
try (axial and radial trapping period = 45 to 210 ms, 
T± = 14 ms) and tunable aspect ratio A = T^/T^ [10]. A 
green laser beam far detuned from the atomic resonance 
is masked to shine on one half of the superfiuid. In a 
time t, the applied potential /7, as experienced by a sin- 
gle fermion, advances the phase of the superfiuid order 
parameter in the exposed region by = 2Ut/h relative 
to the unexposed region. The time t ^ 35 /is is exper- 
imentally adjusted in order to create one high-contrast 
soliton. 

In the strongly interacting regime, the soliton does not 
cause a density depletion within our resolution (see Sup- 
plemental Material). However, it is tied to a phase twist 
in the pair wavefunction. As in the case of vortices [29 , 
the pair wavefunction can be directly observed via a 
rapid ramp to the BEC-side of the Feshbach resonance. 
The ramp converts large fermion pairs into tightly bound 
molecules, empties out the soliton cores and increases the 
soliton width to the final healing length cx 1/ ^Jf^a]^ where 
a/ is the scattering length at the final magnetic field. The 
rapid ramp followed by time of fiight expansion thus en- 
hances the soliton contrast and acts as a magnifying glass 
(for details see Supplemental Material). 

Figure 1 c)-d) reports the observation of solitons in 
a fermionic superfiuid prepared close to the 832 G Fesh- 
bach resonance, at 815 G, for various hold times follow- 
ing the phase imprint. Here, the interaction parameter 
at the cloud center is l/kpa = 0.30(2), where a is the 
scattering length. Figure Ic) shows the optical density 
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FIG. 1. Creation and observation of solitons in a 
fermionic superfluid. (a) Superfiuid pair wavefunction 
A(^) for a stationary soliton, normalised by the bulk pair- 
ing gap Ao, and localised Andreev bound state density ver- 
sus position z, in units of the BCS coherence length ^. (b) 
Schematic of the experiment. A phase-imprinting laser beam 
twists the phase of one half of the trapped superfluid by ap- 
proximately TT. The soliton generally moves at non-zero ve- 
locity Vsoiiton- (c) Optical density and (d) residuals of atom 
clouds at 815 G, imaged via the rapid ramp method [29] , 
showing solitons at various hold times after creation. One pe- 
riod of soliton oscillation is shown. The in-trap aspect ratio 
was A = 6.5(1). (e) Radially integrated residuals as a function 
of time revealing long-lived soliton oscillations. The soliton 
period is Tg — 12(2) T^, much longer than the trapping period 
of Tz — 93.76(5) ms, revealing an extreme enhancement of the 
soliton's relative effective mass M* /M. 



in absorption images taken after time of flight and the 
rapid ramp to ~ 580 G, while Fig. Id) displays residu- 
als obtained by subtracting a smoothened copy of the 
same absorption image. The optical density contrast is 
about 10% (see Supplemental Material). A sequence of 
radially integrated residuals as a function of time is dis- 
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played in Fig. 1 e), demonstrating the soliton to be stable 
for more than 4 s or 100 000 times the microscopic time 
scale h/Ep. This establishes that solitons in fermionic 
superfluids can exist as stable and long-lived excitations. 

The solitons are observed to undergo oscillations in 
the harmonically trapped superfiuid, demonstrating their 
emergent particle nature. However, their period of oscil- 
lation Ts is about one order of magnitude longer than the 
trapping period for single atoms. This directly indi- 
cates an extreme enhancement of their effective mass M* 
relative to their bare mass M = Ngm < 0, where \Ns\ 
is the number of atoms missing in the soliton. Indeed, 
the soliton mass times acceleration M^u^z for periodic 
motion with frequency ujs must be provided by the trap- 
ping force experienced by the missing atoms at position 

NsTtiujIz. This yields a direct relation between the 
relative effective mass M"" /M and the normalised soliton 
period Ts/T^: 



(2) 



In general, the difference between the effective mass M* 
and the bare mass M of the soliton arises from the phase 
slip A(/) across the soliton, which implies a superfiuid 
backflow p!2l. For the soliton to move, an entire sheet 
of atoms thus has to flow past it. The difference M — M* 
is the mass of that sheet, given by the mass density mul- 
tiplied by the entire soliton volume. In contrast, the soli- 
ton's bare mass M is only due to the mass deficit of \Ns\ 
atoms and can become much smaller than M* when the 
soliton is filled. For weakly interacting BECs, where soli- 
tons are devoid of particles, the effective mass is still on 
the same order of the bare mass, (M*/M)gg(^ = 2. This 
leads to an oscillation period that is only y/2 times longer 
than [25l |26l [30l |3l] . In the BCS limit, where only a 
minute fraction Ao/Ep of the gas contributes to Cooper 
pairing, \Ns\ ex Aq/Ef cx e '^f'pi^i and thus the soliton's 
relative effective mass can be expected to become expo- 
nentially large. 

Indeed, as shown in Fig. 2, we find that the soliton 
period, and hence the relative effective mass, increases 
dramatically as the interactions are tuned from the limit 
of Bose-Einstein condensation (Fig. 2a) towards the BCS 
limit. At 700 G, where l/kpa = 2.6(2), the system rep- 
resents a strongly interacting Bose gas of molecules ^ . 
The soliton period is Tg = 4.4(5)T2, already three times 
longer than in the case of a weakly interacting BEC. At 
the Feshbach resonance (Fig. 2d), we measure a soliton 
period of Tg = 14(2) T^, corresponding to a relative effec- 
tive mass of M*/M = 200(50). This is more than fifty 
times larger than the result of mean field Bogoliubov- 
de Gennes theory in three dimensions [121 HO] that pre- 
dicts M*/M = 3. Note that the superfiuid is fully 
three-dimensional: on resonance, the chemical potential 
/i ^ 25fico'x, where uj± is the radial trapping frequency. 
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FIG. 2. Soliton oscillations in the BEC-BCS crossover. 

Shown are soliton oscillations in a trapped fermionic super- 
fluid for various magnetic fields B around the Feshbach reso- 
nance. The soliton period is observed to dramatically increase 
as the system is tuned from the BEC regime (a) to the Fes- 
hbach resonance (d). The measured period, magnetic field 
B and interaction parameter at the cloud center l/kpa were 
(a) Ts/T, = 4.4(5), 700 G, 2.6(2) (b) T./T, = 7.5(9), 760 G, 
1.4(1) (c) Ts/T, = 12(2), 815 G, 0.30(2) (d) Ts/T, = 14(2), 
832 G and 0. The initial atom number per spin state, its decay 
rate and Thomas- Fermi radius range from A^o = 1.1 x 10^, r — 
1.2(2) s and Rty = 135 /xm at B = 700 G to A^o = 2.3 x 10^ 
T = 12(1) s and Rty = 200 /xm on resonance. The aspect 
ratio is A = 6.2(7). Note that dX B — 700 G, the superfiuid is 
short lived due to enhanced three-body loss. At 760 G (b), the 
soliton survived for more than 6 s, comparable to the lifetime 
of the superfiuid at that field. 



Still, for very elongated traps, one expects to reach a 
universal quasi- ID regime where the tight radial confine- 
ment is irrelevant for propagation along the long axis [32] . 
This prompted us to study the dependence of the soliton 
period on the aspect ratio of our trap. 

Figure 3 summarizes our measurements for the soliton 
period and the relative effective mass as a function of 
the interaction parameter l/kpa throughout the BEC- 
BCS crossover, for aspect ratios A = 3.3, 6.2 and 15. 
The strong increase of M*/M towards the BCS regime 
is observed for all trap geometries. The normalised soli- 
ton period Tg/T^ appears to converge to a limiting value 
for the most elongated trap: The normalised period 
changes by only 15% as the aspect ratio is increased by 
more than a factor of two from 6.2 to 15. This indi- 
cates that the soliton dynamics approach the universal 
quasi- ID limit. Even in a much less elongated trap with 
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FIG. 3. Soliton period and effective mass versus in- 
teraction strength in the BEC-BCS crossover. The 

soUton period is shown as a function of the interaction param- 
eter l/kpci in the cloud center, for three different trap aspect 
ratios: A = 15(1) (black circles), 6.2(7) (red diamonds) and 
3.3(1) (orange squares). The soliton period strongly increases 
from the BEC-regime towards the Feshbach resonance, where 
Ts/T, = 12(2) for A = 15(1), and to the BCS side. This di- 
rectly reflects an extreme enhancement of the relative effective 
mass M*/M = T^^/T^, which we attribute to strong quantum 
fluctuations and the filling with Andreev bound states. The 
result for a weakly interacting EEC, Tg/Tz — a/2, is shown as 
the dashed line. The star marks the mean field prediction 
M7M = T'i/T'^ = 3. 



A = 3.3(1) the soliton period is only slightly increased 
by about 30% compared to A = 6.2, accompanied by an 
increased susceptibility of the soliton towards bending or 
"snaking" |9l[T8l|23] (see Supplemental Material). 

We attribute the large relative effective mass M*/M in 
the strongly interacting regime to the filling of the soli- 
ton with uncondensed fermion pairs resulting from strong 
quantum fluctuations. These can also reside inside the 
soliton p!4HT9] , in addition to fermionic Andreev bound 
states. A substantial filling of the soliton will reduce the 
number |A^s| of atoms missing inside the soliton, therefore 
considerably weaken the restoring harmonic force from 
the trap and strongly increase M*/M. Mean- field the- 
ory for the BEC-BCS crossover heavily underestimates 
the role of quantum fluctuations already on the BEC side, 
where it predicts a fraction of uncondensed bosons that 
scales as na^ instead of the correct Vna^ scaling [10]. 
Our experiment thus directly reveals the importance of 
beyond mean field effects for the dynamics of strongly 
interacting fermionic superfluids. 

To demonstrate that the slow soliton oscillations are a 
truly quantum effect and not due to the finite tempera- 
ture of our gas, we investigated the soliton motion as a 
function of temperature for the unitary Fermi gas at the 
Feshbach resonance (see Fig. 4). The soliton period is 



found to be insensitive to changes in temperature within 
the measurement uncertainty (Fig. 4g). The stability of 
solitons is however strongly affected by the thermal frac- 
tion. At low temperatures, the soliton oscillation occurs 
essentially without energy loss, demonstrating dissipa- 
tionless fiow (Fig. 4d and h). For increasing temperature 
we observe anti-damping of soliton oscillations (Fig. 4e 
and h). This is characteristic for a particle with negative 
mass that can lower its energy by accelerating. The en- 
ergy loss is likely due to collisions with thermally induced 
phonons [18^. At even higher temperatures, the soli- 
ton's position becomes less reproducible and its lifetime 
is strongly reduced (Fig. 4f and i). Concurrently, we ob- 
serve increased axial fluctuations in the superfluid (Fig. 
4a-c), some of which appear to have comparable contrast 
to the imprinted soliton. These additional solitons might 
be "thermal solitons" , predicted to occur even in equilib- 
rium in weakly interacting Bose condensates [33]. Simi- 
lar to vortex-anti- vortex pairs in two dimensions, soliton- 
anti-soliton pairs can be expected to spontaneously break 
in ID and proliferate. Note that on resonance, the fastest 
solitons we observe move at the exceedingly slow speed 
of 0.45 mm/s or 5% of the (independently measured) 
speed of sound on resonance. Their sudden disappear- 
ance observed e.g. in Fig. 4f can thus not be related to 
motion close to the Landau critical speed. Instead, their 
decay might be tied to inelastic collisions with thermal 
solitons, as soliton collisions have been found to become 
increasingly inelastic towards the BCS-side in theoreti- 
cal simulations [13 . Another possibility for their decay 
at such low speeds is that the soliton's energy dispersion 
has a minimum at an unexpectedly small fraction of the 
critical velocity [13 . 

In conclusion, we have created and observed long-lived 
solitons in a strongly interacting fermionic superfluid. 
Their period of oscillation and thus their relative effec- 
tive mass dramatically increases as the interactions are 
tuned from the BEC-limit of tightly bound molecules 
towards the BCS limit of long-range Cooper pairs. 
This signals strong beyond mean field effects, likely 
due to uncondensed fermion pairs filling the soliton, 
in addition to purely fermionic Andreev bound states. 
An exciting prospect is to directly study the Andreev 
bound states spectroscopically [TTJ [2j. While they 
are not topologically protected, their lifetime should 
equal that of the soliton, many seconds or 100 000 Fermi 
times, so that they might become a useful quantum 
resource. In the presence of spin imbalance, the soliton 
represents a limiting case of the long-sought Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state of moving 
Cooper pairs [2TJ|22]. Indeed, it is energetically favor- 
able for an excess fermion to reside inside a soliton rather 
than inside the bulk superfiuid. While it is difficult to 
realize the FFLO state in equilibrium, direct engineering 
of soliton trains might produce a long-lived metastable 
analog. 
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FIG. 4. Soliton motion in the unitary Fermi gas at various temperatures, a-c show representative optical densities 
and residuals of the superfluid after the rapid ramp, and d-f the trajectory of solitons for thermal fractions 7(2)%, 9(2)%, 
and 15(3)%, from top to bottom. While at low temperatures, the soliton is the only significant density variation, at higher 
temperatures transverse stripes appear that we interpret as thermal solitons. g shows the soliton period versus thermal fraction, 
demonstrating that the period is insensitive to temperature. The 1/e anti-damping time Ts (h) and the soliton lifetime (i) are 
however strongly dependent on temperature. The soliton lifetime is defined as the time when the probability of observing a 
soliton decreased to 50%. 



Methods Summary 

Preparation. The atomic gas is composed of a bal- 
anced mixture of the two lowest hyperfine states of ^Li 
initially prepared at 760 G [10]. The trapping potential 
results from the combination of a magnetic field curva- 
ture providing harmonic confinement in the axial direc- 
tion and an optical dipole trap (wavelength 1064 nm) 
providing tighter radial confinement. The axial periods 
are =210 ms, 95 ms and 45 ms, respectively, for the 
three aspect ratios considered here. 

Phase imprinting. A step-like intensity profile is 
imprinted on a laser beam (wavelength 532 nm, power 
~200 mW) by means of an opaque mask. A 3 jam reso- 
lution imaging system projects the intensity distribution 
at the mask location onto the atoms (the beam waist at 
the atoms is 60 /im). A phase twist of tt corresponds 
to a pulse time of about 30 /is, much shorter than the 
time scale h/ ja associated to a typical chemical potential, 
on the order of few 100 /is. The pulse duration is finely 
adjusted to yield exactly one soliton with high contrast 
observed after the rapid ramp. Since we imprint only a 
phase step but not a density depletion, sound waves must 
be generated in addition to the soliton [71 [8] . The sound 
waves are found to die out in a quarter axial trapping pe- 



riod when they have reached the edge of the atom cloud. 
For the data above 760 G the soliton is created at 760 G 
and the magnetic field is subsequently ramped (in about 
2 T^) to the final magnetic field where the soliton motion 
is studied. For final magnetic fields below 760 G the soli- 
ton is created at that field. We found that solitons can 
be created directly at the Feshbach resonance as well. 

Thermometry Thermometry from fits of density pro- 
files at the Feshbach resonance to the known equation 
of state [34 yielded an upper limit of T = 0.05T;? = 
10 nK ^ 3Jiw± for our lowest temperatures. At such low 
temperatures far below the critical temperature for su- 
perfluidity, thermometry via fitting is less sensitive to 
small changes in temperature [10]. We therefore use the 
thermal fraction of molecules after a rapid ramp to the 
BEC-side of the Feshbach resonance as a robust ther- 
mometer [To] . 

SUPPLEMENTARY INFORMATION 
Imaging Sohtons 

Solitons gradually fill in as the interaction strength is 
tuned from the BEC-regime to the BCS-regime of the 
crossover. Indeed, in the BCS regime only a minute frac- 
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FIG. 5. Imaging solitons. a Optical density, b integrated ID profiles and c corresponding residuals of a fermionic superfluid, 
prepared at 832 G, after expansion and rapid ramp to various final magnetic fields Bmin- Without any ramp, the superfluid at 
832 G, observed after 9 ms time of flight, does not show a clear signature of the soliton. A 10 ms ramp to the BEC-side before 
expansion at 760 G reduces interactions but still only reveals a very faint trace of the soliton. For Bmin < 700 G the soliton is 
revealed, d Sequence of the magnetic field ramp, indicating time of flight (TOF), final ramp field Bmin and the imaging pulse 
at 760 G. e The maximum depletion detected in the optical density, in units of the standard deviation a found outside the 
soliton. The detection threshold of 2.5 a is indicated. 



tion Aq/ Ef of the gas is Cooper paired, and only this 
fraction is missing at the soliton's center, where the pair 
wavefunction A(z) is maximally depleted. The contrast 
in the particle density thus vanishes. Indeed, absorp- 
tion images of expanded atom clouds at the Feshbach 
resonance (Fig. 5a) do not show any observable 3%) 
contrast. However, the modulus of the pair wavefunction 
itself can be imaged by a rapid ramp technique similar 
to what was used for the observation of vortex lattices 
in the BEC-BCS crossover [l^ A magnetic field 

ramp to the weakly interacting BEC-regime turns large 
fermion pairs into molecules. An absorption image of 
molecules thus approximately reflects the magnitude of 
the fermion pair wavefunction before the ramp. In addi- 
tion, the ramp reduces the interaction strength and thus 
increases the coherence (healing) length of the superfluid, 
which increases the soliton contrast and increases their 
width. 

The rapid ramp is illustrated in Fig. 5. Starting for ex- 
ample with a superfluid at the Feshbach resonance, the 
magnetic field is first quickly ramped over 10 ms to 760 
G, on the BEC-side of the Feshbach resonance where in- 
teractions are weaker and fermion pairs are more tightly 
bound. Next, the cloud is released from the trap. After 
1 ms, the magnetic field is rapidly ramped over 1 ms 



to ~ 580 G, where interactions are essentially absent 
and fermion pairs have fully turned into tightly bound 
molecules. The molecular cloud further expands for 4 ms 
at 580 G, after which the magnetic field is re-ramped over 
1 ms to 760 G. After an additional 2 ms of expansion at 
760 G, the molecules are imaged via absorption imaging. 

Solitons can be identified easily in the absorption im- 
ages by eye. However, to automize soliton detection we 
implemented the following method. For each absorption 
image, we first generate a residual profile by subtracting 
a smoothened version of the optical density profile from 
the actual optical density profile. We determine the stan- 
dard deviation of fluctuations a and identify a depletion 
in the residual profile as a soliton if its depth is greater 
than 2.5 a. 

We have found the rapid ramp technique necessary to 
reveal solitons in the strongly interacting regime, which 
is another indication of their strong filling, next to their 
slow period and enhanced relative effective mass. To 
show the importance of the rapid ramp, we have varied 
the final field of the rapid ramp ^min between 500 and 
832 G. The depth of the maximum depletion, normalized 
by cr, is shown in Fig. 5. For ramp fields 5min < 650 G 
solitons are clearly revealed. 



FIG. 6. Observation of the snake instability in a fermionic superfluid. a-c Absorption images exhibiting "snaking" 
sohtons in expanding fermionic superfluids. d-f Same as a-c) but with a white dashed hne as guide to the eye. The superfluids 
were prepared at B = 760 G in a trap with aspect ratio of 2.9 (a) and 3.1 (b-c) with an axial trapping period of Tz = 44 ms. 
The images were taken 200 ms after the phase imprint of the soliton, using the standard rapid ramp to Bmin = 580 G and a 
total time of flight of 9 ms. 



Snake Instability Planar solitons in three dimensions 
are not only thermodynamically unstable towards accel- 
erating, but also dynamically unstable towards shape ex- 
citations - the so-called "snake" instability [9l [181 ES]. 
This excitation of the soliton plane along the radial di- 
rection grows exponentially until the soliton decays into 
vortices. This is suppressed in elongated trap geometries. 
In weakly interacting BECs, solitons are expected to be- 
come dynamically unstable when the chemical potential 
(1 of the condensate becomes larger than « 2AhjO^^ where 
uj±_ is the radial trapping frequency [35]. The solitons in 
our strongly interacting fermionic superfluid appear to 
be much more robust, as they are long-lived even when 
the chemical potential ji ~ 2bhx)^. Still, we were able to 
observe the snake instability by reducing the trap aspect 
ratio to about 3. Examples are shown in Fig. 6 where the 
depletion revealed by the rapid ramp no longer follows a 
straight line but rather a wavy trajectory, characteristic 
of the snake instability p [T8ll23] . 
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